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Abstract 

We consider the totally asymmetric exclusion process (TASEP) in one 
dimension in its maximal current phase. We show, by an exact calculation, 
that the non-Gaussian part of the fluctuations of density can be described 
in terms of the statistical properties of a Brownian excursion. Numerical 
simulations indicate that the description in terms of a Brownian excursion 
remains valid for more general one dimensional driven systems in their max- 
imal current phase. 
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1 Introduction 



Exclusion processes |2J El with open boundaries have attracted much 
attention as simple models of an open non-equilibrium system in contact 
with two reservoirs having different chemical potentials [SJ |S] • Despite 
their simplicity, these models exhibit properties believed to be characteristic 
of realistic non-equilibrium systems, such as long range correlations [7JIB1IH1 
1101 lllj and phase transitions in one dimension [IJ [5J • 

A number of exact results have been obtained for the one dimensional 
exclusion process with open boundaries, using the fact that the weights of the 
microscopic configurations in the stationary state can be calculated exactly 
0El[I2[I3imiI3l]iIll]I!- The goal of the present paper is to show that 
the density fluctuations of the totally asymmetric exclusion process in its 
stationary maximal current phase can be expressed in terms the statistical 
properties of a Brownian excursion. 

A Brownian excursion ^Sj is a Brownian path Y{x) conditioned to re- 
main positive, i.e. such that Y{x) > for < x < 1, with the boundary 
condition 1^(0) = 1^(1) = 0. In the following we will consider the Brown- 
ian excursion Y{x) normalized in such a way that the probability density 
of being at heights yi,y2,--Uk at times < x\ < %2 < ---Xk < 1 (i.e. that 
Y{x\) = j/i, ...Y(x k ) = y k ) is given by 

^ fc-i 
Prob(2/i,...yfc;aJi,...a; fc ) = —= h(y X ]Xi) h(y k -\-x k ) TT g(y p ,y p+1 ;x p+1 -x p ), 

;=i 

(i.i) 

where the functions h and g are defined by 



and 

g(y,y';x) 



Hy;x) = -^j2 e " (I- 2 ) 



(1.3) 



One can check that is normalized, i.e. that 

f'OO roc 

/ dyi... dy k Pvob{yx,...y k ;xi,...x k ) = 1, 
Jo Jo 

using the identities 

1 f 00 

— dy h(y; x) h(y; 1 - x) = 1, 

dy' /i(y'; x) y; a/) = h(y; x + x'), 

dy" g{y, y";x) g(y" y'; x') = g(y, y'; x + x'). 

The one dimensional totally asymmetric simple exclusion process (TASEP) 
with open boundary conditions is defined as follows: each site i (with 
1 < i < L) of a one dimensional lattice of L sites is either occupied by 
a single particle or empty, and the system evolves according to the following 
continuous time dynamics: if a particle is present on site i (for 1 < i < L— 1), 
it hops at rate 1 to site % + 1 if this site is empty. At the left boundary, site 
1 is filled at rate a by a particle if it is empty. At the right boundary, if a 
particle is present at site L, it is removed at rate (3. Each microscopic con- 
figuration can be described by a set of L binary variables Tj, the occupation 
numbers (t^ = 1 if site i is occupied and r» = if it is empty). When a 
and f3 lie in the interval (0, 1), the case we shall be concerned with here, the 
input and exit rates a and (5 can be thought as resulting from the system 
being in contact with a left and and a right reservoir at densities a and 1 — /3 
respectively [TT] . 

In the steady state, which is unique for such a Markov process, one can 
try to determine correlation functions, which we will denote (r^r^.-.r^). 
One can also divide the L sites into k boxes of L\, L2, ■■■L k sites and try to 
determine the probability that in the steady state N\ particles are present 
in the first box, A^2 in the second box,... N k in the fcth box. 
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In the present paper we are going to show that in the maximal current 
phase @1 El of the stationary TASEP, corresponding to 

a > i and > ^ (1.4) 

the correlation functions of the occupation numbers Tj in the steady state 
are given for large L by 

> - 1) ( r - - 1) ••■ ( r - - £)) * *n ^dxt..^ 



(1.5) 

where %\ < 12, ■ ■■ < i k and the right hand side is to be evaluated at x p = i p /L. 
The averages are taken with respect to Ijl.ljl . 

Our second result concerns the fluctuations of the numbers Ni,...N k 
of particles in boxes of length L\ = L{x\ — Xo),L,2 = L(x2 — x{), ...L k = 
L(xk — Xk-i) with xo = < x\ < .... < Xk-i < Xk = 1. Define fi p to be the 
rescaled fluctuations of the number of particles, in box p, 

^ = N p -L{ Xp -_x p ^)/2 (16) 

We are going to show that their probability density Q(fii, ■■■fJ-k) • •■%k-i) 
is given for large L by 

Q{m, ...fik;xi, ...acfc-i) = / dyi... / dy k -i>< (1.7) 

-(tyi + yi-i -Vi) 



k „ r /«. . x 2 



Prob(yi, xi, ...x fc _i) TT - == exp 



where y = y k = 0, x = 0, x fc = 1 

As we shall see in section 2 the product in the integrand of (|1.7|) is just 
the conditional probability of /ij, /z& given the values of the 1" process, 
Y(xi) = yi, ...y(xfc_i) = 2/fc-i- This conditional probability is just a product 
of Gaussians with means (yi — j/j_i)/2 and variances (xj — Xj_i)/8. Since 
(|1.7jl is valid for arbitrary number and sizes of boxes it is equivalent to 
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the statement that the "fluctuation field" of the particle density p(x) in 
the maximum current phase can be written as a sum of two independent 
processes, 



VI 



p(x) 



B(x)+Y(x) 



Here p(x) is the empirical density at x, defined by 



VI 



P /p{x) -\) dx ^ VI 



7f E 



i=x p — iL 



while Y(x) is the (generalized) derivative of the Brownian excursion de- 
scribed by B{x) is a white noise, the derivative of a Brownian path, 
normalized so that 

' " (1-9) 



\B[x) - B{x')] 2 



\x — X 



and B and Y are independent. 

For the integrated fluctuation r(x,x') of the density in the macroscopic 
segment (x, a/), 



r(x, x') 
one has 



, . , , l\ B(x') + Y(x') - B(x) - Y(x) 
' { P{V) ~ 2 J " 2VI ( } 



r(x, x') 



r\2 



[B{x) - B{x')} 2 + [Y{x) - Y{x')} 2 



4L 



;i.ii) 



If one considers now the fluctuation in a very small segment away from the 
end points, then Y(x') — Y{x) behaves like B(x') — B(x) and so 



r(x, x') 



A2 



\X — X\ 

4L 



for \x — x'\ — > 0, 



indicating that locally the measure is Bernoulli. On the other hand if one 
considers the particle number fluctuation in the whole system 



r(0,l) 2 ~ 



1 

8l' 
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which means that the fluctuation of the total number of particles is one half 
of what it would be for a Bernoulli measure at density 1/2 

One can check that (|1,5|) extends to arbitrary correlations what was 
already known for the one-point and two-point functions (equations (47) 
and (52) of [S]) when a = (3 = 1. Also (|1.7|) extends to an arbitrary number 
of boxes the result (6.15) of jllj valid for a single box. 

Our derivation presented in section 2 is a generalization of the method 
used in Numerical simulations reported in section 3 indicate that the 
description in terms of a Brownian excursion remains valid in the maxi- 
mal current phase of other models for which the steady state is not known 
exactly. 

2 Derivation of flPft and ( IT77I) 

2.1 The matrix method 

For the steady state of the open TASEP described in Section ^ the proba- 
bility P{{t{\) of any microscopic configuration {t{\ can be written as ^2] 

P({Tj}) _ M 

where the normalization factor Zl is given by 

Z L = (W\(D + E) L \V), (2.2) 

and the matrices D, E and the vectors {W\, \V) satisfy the relations, 

DE = D + E, (2.3) 
/3D\V) = \V), (2.4) 
(W\aE = (W\ . (2.5) 
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Prom the algebra (|2.1|) ~ ()2.5j) all equal time steady state properties can 
be calculated. For example the average occupation (t$) of site % is given by 

(n) = y ' ( 2 - 6 ) 

and the two point function is, for i < j, 

(W\(D + Ef- 1 D(D + Ey- 1 - 1 D(D + E) L ^\V) 
\ T i T j) = — • (2.7) 

The probability of finding N%, .JVj. particles in subsystems of lengths 
L\, ...Lk can be written as 

(N N) (WlX^m X L2 (N 2 )... X Lk (N k )\V) 

quM-Lkl^i, -Nk) = -g- (2.») 

where Xl(N) is the sum over all the products of L matrices containing 
exactly N matrices D and L — N matrices E. This can be written as 

X L (N)=J d6 (De 2 ™ e + EJ e ~^e ^ ^ 

The algebraic rules (|2.3H2.5j) allow one to calculate all the matrix elements 
appearing in ()2.1I2.2I2.6H2.8|) without using any explicit representation of 
the matrices D and E or of the vectors (W\ and \V). Working with a 
particular representation might be convenient but of course the steady state 
properties, such as correlation functions and current, do not depend on the 
particular representation used. 

To derive the expressions l)1.5|1.7|) we find it convenient to use a par- 
ticular representation of 1)2. 312. 5 Jl (which was already used in section 6.3 of 
d): 

00 

D = ^2\n)(n\ + \n)(n+l\, (2.10) 

n=l 

00 

E = ^ \n){n\ + \n + l)(n|, (2.11) 

n=l 



where the vectors |1), |2), ...|n)... form an orthonormal basis of an infinite 
dimensional space (with (n\m) = 5 ntm ). In this basis, the vectors \V) and 
(W\ satisfying (|2.4j) . l|2.5jl are given by 



oo 



\V) = ±^[ L ^-) n \n), (2-12) 



n=l 



°° / 1 — \ n 

W=E — H (2-13) 



, a 

n=l 



In the following we will assume that 

l<a + (3 , a < 1, and < 1. (2.14) 

so that (W|V) and all the matrix elements (W|Xi-X2...|V) are finite and 
positive when the matrices X\,X2- are polynomials of matrices D and E 
with positive coefficients. This condition is not the same as the condition 
1|1-4J) of being in the maximal current phase. We will see later (|2.26f) how 
condition Q1.4JI appears. 

Note that as long as L is finite, all the matrix elements are rational 
functions of a and (5 and so all the expressions derived assuming (|2.14f) 
could be analytically continued to the whole range of values of a and (3. 

2.2 The sum over walks 

Let us introduce the set of discrete walks w of L steps which, at each step, 
either increase by one unit, decrease by one unit, or stay constant, with the 
constraint that they remain positive. Each such walk w can be described by 
a sequence of L + 1 integers {n,(u>)} satisfying for all i = 0, 1, L 

rii > and \rii — < 1 

To each such walk u>, one associates a weight tl(w) defined by 



where v(n, n') is given by 



'2 if |n - n'\ = 
v(n, n') = < 1 if \n — n'\ = 1 
,0 if |n — n'| > 1 

It is easy to check from (j2.10l2.11j) that for n > 1 and n' > 1, one has 

v (n, n) = {n\D + E\n) 

and it follows that 



(W\(D + E) L \V)=J2^)- 

w 

These weights define a measure v{w) on the walks w 



(2.16) 



v{w) 



(2.17) 



It follows from (12.1 012.1 1 ft that 
(n\D\n'} 



(l + n' -n)(n\D + E\n') 



which combined with ()2.7|) yields, for i < j, 

( T i T j) = 7 X] t 1 + n ^ W ) ~ "i-lWl I 1 + n j( W ) - n j-l( W )] • ( 2 - 18 ) 

ui 

More generally, for i\ < %2 < ... < ik, 

(n 1 --n k ) = 2kJ2^(w)[l + n h -n il - 1 ]....[l + n ik - n ifc _i] , (2.19) 

(where to avoid heavy notation we have not repeated the w dependence of 
all the n, 's). This can be rewritten as 



1 



(2.20) 



which is the exact finite L version of Q1.5|) . 

The expressions of (n\(D + E) L \ri) and of {n\X L {M)\ri) defined in (123]) 
are known (see equations (6.24) and (6.65) of which had been derived 
by recursions over L in |2*U1 . 

l, m (2L)\ (2L)\ 



(n\(D+E) L \n') 



(L + n- n')\ (L + ri - n)\ (L + n + n')\ (L - ri - n)!' 

(2.21) 



(n\X L (N)\n') - 



(N)\ (L - N)\ (N + n - n')\ {L-N -n + n')\ 
(LI) 2 



(N + n)\{L-N - n)\ (N - n')\ (L - N + n')V 

where any negative factorial is defined to be infinite (i.e. the matrix elements 
are non-zero only when N — L < n' — n < N < L). 
Let 

(n\X L {N)\n') 

FL,N(n,n)- {nl{D + E)Lln ,y (2-23) 

The probability qL 1 ,L 2 ...L k (Xi, ..N^) that there are Ni particles in the first 
L\ sites, N% in the next L2 sites, etc., is then given by 

k 

qL u L 2 ...L k (Nt,..N k ) = ^v(w)l[F Lt , Ni (n Mi _ 1 (w),n Mi (w)) (2.24) 

w i=l 

where Mq = and Mi = Mj_i + Li. This is the exact finite L version 
of (|1.7|) . It shows clearly that given w the {Ni} are independent random 
variables. 

2.3 Derivation of flUD 

Let us first evaluate for large L the normalization factor ()2.16j) 

00 00 /-, \ n /l R\ n ' 

E n W = EE( — J (V) W^ + ^V) (2.25) 

«) n=\n'=\ ^ / \ P / 
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Using Q2.21j) (or its large L behavior easily obtained by the Stirling formula) 
one can show that if 

and /3 > -, (2.26) 

then the walks which dominate the sum for large L are those which have 
both no(w) and til(w) of order 1 and which remain at distances of order 
L 1 / 2 from the origin. Condition Q2.26JI corresponds to the maximal current 
phase |121 119j . It is more restrictive than Q2.14J1 which assures only that 
(W|V) is finite. In the range where Q2.14JI is satisfied but (|2.26|) is not, the 
walks which dominate (|2.25|) are walks such that either no(w) or til(w) is 
of order L, and where the Brownian excursion picture does not apply. 

The large L expressions of matrix elements of the form ()2.21|) can be 
easily obtained using Stirling formula and one gets, if and n and n' are of 
order s/L, 



(n\(D+E) L \n') ~ 



4 L {-(n-n') 2 \ (-(n + n 1 ) 2 



exp — exp 



L / "V L 



. (2.27) 



ML 
Hence for 

i p = Lx p , n ip = L 1/2 y p , (2.28) 

one gets 

4 l ( x p+i- x p) 

{n ip \(D + Eyp+i % *\n ip+1 ) ~ = g{y p , y p+1 ; x p+1 - x p ) 

V L 

where g is defined as in . This formula remains valid even if n j and / or 
n i P +i are OI " order 1. For example one obtains that way that if rii p is of order 
1 and n ip+1 ~ \ZZy p+1 , 

^L(x p+ i-Xp) 

(n ip \{D + E) l ^-^\n lp+1 ) ~ — 2n ip h{y p+1 ; x p+l - x p ) 

so that (|2~^oT) and(j22Z|) give 

(WUD + E) L \V) ~ 4L+1 {1 ~ a)a <± ~ W (2 29) 
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The correlation function n^.-.n^ of the heights of a walk w at positions 
il, ...ik is then given by 



n ,n il ,...n ik ,n L 

(no|( J D + E)«K)....K|p + E) i -Hn i )(^) no (i/ 



and for i\ = Lx\, = Lxfc one gets 

n— n~ = L k l 2 Y(x 1 )...Y(x k ) 
which using (f2~2~U)) leads to (|T3)) . 

2.4 Derivation of (TT7fI> 

For large L, with n, n' and N — L/2 of order \/X> one can easily see from 
(t2~2~2l that 



(n\X L (N)\n) ~ 2 



ii 

vrL 



qA±±±v/_ o V 

e £ £ 



e z l l 



where AA^ = N — L/2, which can be rewritten as 



— (2AiV+n — n ) 



(n\X L (N)\ri) ~ 2 — e r 

7TL 



e l - e l 



Thus Fi,^(n,n') defined in (|2,23|) becomes for n,n' and AA" of order \fh 



F L ,N(n,n') 



exp 



-(2 AN + n-n 
L 



r\2i 



and this shows that in the large L limit (j2.24j) reduces to fi .7j) . 
2.5 Origin of (TOil 

As already noted in section 1, expression (|1.8jl is essentially equivalent to 
Q1-7J) so the derivation of (|1-7|) above also gives (|1.8jl . It is interesting, 
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however, to understand the origin of Q1.8|) directly from the microscopic 
picture involving the walks w. To do that let us define the joint distribution 
u(w,t) of w and r = {t{\. It follows directly from (|2.19|) that 

v(w,t) =v(w)P(t\w), (2.30) 

where v{w) is defined in (|2.17|) and the conditional probability of the {u} 
given w is a product measure 

P( T \w) = J| V 1A '-. (2.31) 

This leads to a simple way of generating steady state configurations of 
the occupation numbers {t{}. First one generates a random walk w of L 
steps according to the measure v{w). Then according to (|2.31|) a steady state 
configuration {rj} is obtained by taking t,- l = 1 whenever rii(w) — nj_i(u>) = 
1, Tj = whenever rii(w) — nj_i(u>) = —1 and by choosing Tj = or 1 with 
equal probabilities for each i such that rii(w) — nj_i(ui) = 0. Therefore the 
fluctuations of the density have two contributions: the random choice of the 
walk w which is at the origin of Y(x) in ()1.8|) and once w is chosen, the 
random choices of the Tj for the flat parts of the walk, which correspond to 
B(x) in (J1.8D . The normalization of B(x) in (|1.9|) arises from the fact that, 
for large L, rii(w) = rii_i{w) for approximately half the steps of the walk w. 

3 Simulations 

An interesting question is to know whether the fluctuations of density of 
one dimensional driven diffusive systems in their maximal current phase 
take always the form (jl.8|) . once properly normalized. 

First, we believe that our results also hold for the general ASEP where 
particles can also jump to the left with rate q < 1. In this case also, one 
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can choose a representation of the matrices, used recently to calculate the 
large deviation function in the weak asymmetry limit |22| . such that matrix 
elements can be thought as sums over weighted walks which do not cross 
the origin. We will not discuss this further here. 

In order to test whether the results obtained here for the TASEP remain 
valid for a broader class of models, we consider in this section a generalisation 
of the TASEP |23l I24j in which particles are extended. In this model, each 
particle occupies d consecutive sites of a lattice of L + d — 1 sites and the 
exclusion rule forbids that any site of the lattice is occupied by more than 
one particle. There are thus L possible positions for a single particle of size 
d on the lattice of L + d — 1 sites. By convention, we say that a particle 
is at site % when it covers sites i,i + 1, ...i + d — 1. The system evolves 
according to the following rule: during each infinitesimal time interval dt, 
each particle jumps one step to its right with probability dt provided that 
this is allowed by the exclusion rule. Moreover if the first d sites are empty, 
a new particle is injected at site 1 with probability adt and if a particle 
covers sites L, L + 1, ...L + d — 1, it is removed with probability [3dt. For 
d = 1 the problem reduces to the TASEP discussed in sections 1 and 2. 

For general d, the current Jd and the density pd in the maximal current 
phase are given by [2S] 

J « = W^r (31) 

Pd= Vd(V3 + i) (3 - 2) 

These expressions can be understood by considering a system of M such 
particles on a ring of L sites, and by using the fact that all allowed configu- 
rations are equally likely (see the appendix). 

In this ring geometry, the number m of particles on I consecutive sites 
fluctuates in the steady state and it is possible to show (see the appendix) 
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that for d <C I <C L, one has 

(m 2 ) — (m) 2 



A d (p) (3-3) 



/ 

with 

A d (p)=p(l-p)-p 2 (d-l)(2-dp), 

This expression can also be recovered from the pressure |25| of this system 
p = log ^ - pr^f^- J using Ad(p) = p^- In our simulations we tried to see 
whether, for large system sizes with open boundary conditions, the fluctu- 
ations of density in the maximal current phase would still be given by the 
statistics of a Brownian excursion with (|1.8j) and ([l.lOjl being replaced by 



p(x) -P d = \j [B{x) + (3-4) 

and 

r(x, x') = f X dy (p(y) - p d ) = ^M[B(x') + Y(x') - B(x) - Y{x)\ 

(3.5) 

If we look at the statistical properties of the normalized number s{x) of 
particles in a box of size Lx centered at the middle of the system 



°W = \hA^r(^,±±Z) (3.6) 



we get, using ([3.5)1 and the properties of Y(x 

2x-3x 2 3-2x + 3x 2 _ x (l-x\ 3(1 - x)y/x 



s(x) 2 = + ■ cos" 1 - -i (3.7) 



and 



3x 2 (x - 2)(5x - 2) 15-3x 2 (x-2)(5x-2) _ x (\-x 



s(x) 4 = 1 cos , 

4 47T V 1 + x 

+ \~ X , ^(15x 3 - llx 2 - 25x - 15) (3. 
27r(l + x) 
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Figure 1: The moments s(x) 2 and s(x) 4 of the fluctuations of the number 
of particles s(x) (normalized as in 1)3. 9|) ) between sites and L^-^- for a 

system of L = 601 sites and for 3 sizes of particles (d=l, 3, 5). The curves 
are the analytic predictions (|3.7I) and (|3.8|) . 
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We have simulated this model for three sizes of particles d = 1, 3 and 5 
for an open system of L = 601 sites when a = (3 = 1. Let N x be the number 
of particles in the interval < i < £ "Mp (we count as iV x the number 

of sites in the interval covered by a particle, divided by the length d of one 
particle). The normalized number of particles s(x) is related to N x by 

N x -Lxp d 

six) = . (3.9) 
s/LA^Pd) 

We have measured the second and fourth moments of s(x) in the steady 
state, averaged over typically 10 s updates per site. In figure ^ we compare 
the results of our simulations with the theoretical predictions (|3.7|) and Q3.8[) . 
The fact that the curves for the different choices of d coincide indicate that 
the fluctuations, once properly normalised, are universal (i.e. do not depend 
on d) 



4 Conclusion 

The main result of the present paper is that the steady state fluctuations 
of the density profile in the TASEP can be written in terms of a Brownian 
excursion as in (|1.5J1.7I1.8|) . 

Our simulations of a more general model in which particles are extended 
indicates that the fluctuations of the density profile may be universal. It 
would of course be nice to make other numerical tests of this universality 
and to see whether it could be understood by a more macroscopic approach, 
such as in ISHIEZI- 

Another interesting question would be to see whether the time dependent 
fluctuations of the density profile would arise from some simple stochastic 
dynamics of the Brownian excursion. 
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Appendix 



Let us consider the TASEP on a ring of L sites with A particles of size d 
as in section 3. This hard rod problem has a long history starting with the 
works of Lee and Yang All allowed configurations (i.e. configurations 
which satisfy the exclusion rule) are equally likely and let us call Zl(N) 
their number. 

The number zl{N) of configurations of A particles on a lattice of L 
sites with open boundary conditions is given by, Zq(N) = z\(N) = 22(A) = 
...Zd-i(N) = 5jv,0) an d for L > Nd 

z L (N) = (L " iVd + Ar)! 
L{ ' N\(L-Nd)\ 

(one can easily check that zl(N) satisfies the recursion Zl{N) = zl-i(N) + 

z L - d (N-l)). 

By considering that site 1 on a ring is either empty or covered by one 
particle, one can express Zl(N) in terms of the partition functions of the 
open systems 

The current J on the ring is given by 

z L _ d _!(A-l) _ N(L-Nd) 
Z L (N) L(L — Nd + A — 1) 

For large L and A, at fixed density p = N/L, J becomes 

j_ P(l-pd) 
l-p(d-l) 

which gives the expressions (|3.'2(l and Q3.1JI when the current is maximal. 

One can calculate the correlation function (rjTj + fc) between the occupa- 
tions on the ring 

1 N 

{TiT i+k ) = z k-d{n) ^L-fc-d(A - n - 2) 



Zl{N) „ 

n=0 
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Using the above expressions of Zl(N) and zl(N), one gets for large L and 
N, keeping p = N/L, k and n fixed 

{nn+u) = Y, z k-d(n) p n+2 (1 - pd) k - nd - d [1 - p(d - 

n>0 

then as 

„JV<2 



L>7Vd 

one gets that (for arbitrary e) 



X' 



(1 - x) N +l 

L>Nd v 1 



..2„-de 



-ke P e 

fc>l 



" (1 - pd)(l - e" e ) + p{l - e~ d€ ) 
and as 



~2„-e 

fc>l 



1-e" 



one finds by expanding the above expressions in powers of e that 
S^lrr \ ItMt \ {d-l)p*{2-pd) 

fe>l 

so that the fluctuations of the number m of particles on I consecutive sites, 
with d <C I -C L is given by 

(m 2 ) — (m) 2 



(p-^)-(d- 1)^(2-^ 



as given in 
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Figure Captions 



Figure ^ The moments s(x) 2 and s(x) 4 of the fluctuations of the number 
of particles s(x) (normalized as in (|3.9jO between sites and L for a 

system of L = 601 sites and for 3 sizes of particles (d=l, 3, 5). The curves 
are the analytic predictions (|3.7|) and (|3.8|) . 



20 



References 

[1] H. Spohn, Large Scale Dynamics of Interacting Particles (Springer- 
Verlag, Berlin 1991) 

[2] B. Schmittman, R. K. P. Zia, Statistical mechanics of driven diffusive 
systems (Academic Press, London, 1995). 

[3] T. M. Liggett, Stochastic interacting systems: contact, voter, and ex- 
clusion processes (Springer- Verlag, New York, 1999). 

[4] J. Krug, Boundary-induced phase transitions in driven diffusive sys- 
tems, Phys. Rev. Lett. 67, 1882-1885 (1991). 

[5] B. Derrida, E. Domany, D. Mukamel, An exact solution of a one di- 
mensional asymmetric exclusion model with open boundaries, J. Stat. 
Phys. 69, 667-687 (1992). 

[6] G. Schiitz, E. Domany, Phase transitions in an exactly soluble one- 
dimensional exclusion process, J. Stat. Phys. 72, 277-296 (1993). 

[7] H. Spohn, Long range correlations for stochastic lattice gases in a non- 
equilibrium steady state, J. Phys A. 16, 4275-4291 (1983). 

[8] B. Derrida, M.R. Evans, Exact correlation functions in an asymmetric 
exclusion model with open boundaries, J. Phys. I France 3, 311-322 
(1993). 

[9] B. Derrida, J. L. Lebowitz, E. R. Speer, Free energy functional for 
nonequilibrium systems: an exactly solvable case, Phys. Rev. Lett. 87, 
150601(2001). 



21 



[10] B. Derrida, J. L. Lebowitz, E. R. Speer, Exact free energy functional for 
a driven diffusive open stationary nonequilibrium system, Phys. Rev. 
Lett. 89 0300601 (2002) 

[11] B. Derrida, J. L. Lebowitz, E. R. Speer, Exact large deviation functional 
of a stationary open driven diffusive system: the asymmetric exclusion 
process, J. Stat. Phys. 110 775-809 (2003) 

[12] B. Derrida, M. R. Evans, V. Hakim, V. Pasquier, Exact solution of a ID 
asymmetric exclusion model using a matrix formulation, J. Phys. A 26, 
1493-1517 (1993). 

[13] S. Sandow, Partial asymmetric exclusion process with open boundaries, 
Phys. Rev. E 50, 2660-2667 (1994). 

[14] K. Maffick, S. Sandow, Finite dimensional representations of the 
quadratic algebra: applications to the exclusion process, J. Phys. A 30 
4513-4526 (1997). 

[15] F.H.L. Essler and V. Rittenberg, Representations of the quadratic 
algebra and partially asymmetric diffusion with open boundaries 
J. Phys. A 29, 3375-3408 (1996) 

[16] T. Sasamoto, One dimensional partially asymmetric simple exclu- 
sion process with open boundaries: Orthogonal polynomials approach, 
J. Phys. A 32, 7109-7131 (1999). 

[17] R. A. Blythe, M. R. Evans, F. Colaiori, F. H. L. Essler, Exact solution 
of a partially asymmetric exclusion model using a deformed oscillator 
algebra J. Phys. A 33, 2313-2332 (2000). 

[18] J. Pitman, Combinatorial Stochastic Processes, Lecture Notes from 
Saint Flour Course, July 2002 



22 



[19] V. Popkov, G.M. Schiitz, Steady state selection in driven diffusive sys- 
tems with open boundaries, Europhys. Lett. 48, 257-263 (1999). 

[20] B. Derrida, M.R. Evans, D. Mukamel, Exact Diffusion Constant for 
One-Dimensional Exclusion Models J. Phys. A 26, 4911-4918 (1993) 

[21] K. Mallick, Shocks in the asymmetry exclusion model with an impurity, 
J. Phys. A 29, 5375-5386 (1996). 

[22] C. Enaud, B. Derrida, Large deviation functional of the weakly asym- 
metric exclusion process, in preparation 

[23] G. W. Lakatos and T. Chou, Totally asymmetric exclusion processes 
with particles of arbitrary size, J. Phys. A 36, 2027-2041 (2003). 

[24] L. B. Shaw, R.K.P. Zia, K. H. Lee, Modeling, simulations, and anal- 
yses of protein synthesis: Driven lattice gas with extended objects 
cond-ma t/0302128 

[25] T.D. Lee, C.N. Yang, Statistical theory of equations of state and phase 
transitions II. Lattice gas and Ising model, Phys. Rev. 87, 410-419 
(1952) 

[26] L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, C. Landim, Fluc- 
tuations in stationary non equilibrium states of irreversible processes, 
Phys. Rev. Lett. 87, 040601 (2001). 

[27] L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, C. Landim, 
Macroscopic fluctuation theory for stationary non equilibrium states, 
J. Stat. Phys. 107, 635-675 (2002). 



23 



